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QUANTUM STOCHASTIC CALCULUS AND QUANTUM 
NONLINEAR FILTERING 

V. P. BELAVKIN 


Abstract. A ^-algebraic indefinite structure of quantum stochastic (QS) cal¬ 
culus is introduced and a continuity property of generalized nonadapted QS 
integrals is proved under the natural integrability conditions in an infinitely 
dimensional nuclear space. The class of nondemolition output QS processes in 
quantum open systems is characterized in terms of the QS calculus, and the 
problem of QS nonlinear filtering with respect to nondemolition continuous 
measurments is investigated. The stochastic calculus of a posteriori condi¬ 
tional expectations in quantum observed systems is developed and a general 
quantum filtering stochastic equation for a QS process is derived. An applica¬ 
tion to the description of the spontaneous collapse of the quantum spin under 
continuous observation is given. 


Introduction 

The problem of description of continuous observation and filtering in quantum 
dynamical systems can be effectively solved in the framework of quantum stochastic 
(QS) calculus of nondemolition input-output processes first developed for quantum 
unitary Markovian evolutions in ^ . 

In contrast to classical probability theory, the conditional expectations defining 
a posteriori states of a quantum system with respect to a subalgebra may not exist 
in general and the existence depends on the algebra of observables and on the a 
priori state. 

During the preparation of the measurement of a Quantum System the change 
necessary to produce the a priori compatible state as a mixture of the a posteriori 
states is referred to in quantum physics as the demolition of the system. The 
latter involves a change in the initial state by the reduction of the algebra of the 
system during such a preparation. The nondemolition principle provides a sufficient 
condition for the algebra of the quantum system to be prepared for the measurement 
in an initial state. 

The mathematical formulation of the nondemolition principle for the observ¬ 
ability of a class of quantum processes was given in [2] and investigated in subse¬ 
quent papers Hi]. This fundamental principle of quantum measurement theory 
means that if a QS process Xt is indirectly observable by the measurement of an¬ 
other process U then X and Y must satisfy the one sided commutativity condition 
[Xt,Ys\ = XfYs — YgXt = 0 for alH > s but not for t < s. 
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In the physical language this means that the measurements of Y in real time do 
not demolish the quantum system X (which has been prepared for the observation) 
at the present time or in the future. The condition given above, however, shows 
that, though the past of X (priori to t) can never be observed, it is demolished 
by the observation of the V process. Mathematically it can be expressed as the 
decomposability of the algebra At generated by {-^(s) ■ s > t}, describing the 
present and future of the system with respect to the spectral resolution of any 
Hermitian operator of the algebra B* generated hy {Y{s) : s < t}. 

In this paper we show using method of quantum filtering that the nondemolition 
condition given above is necessary and sufficient for the evaluation of a posteriori 
mean values oi X G At given are arbitrary initial state. In other words we prove 
that a quantum system is statistically predictable by a measurement procedure, iff 
the observable process satisfies the nondemolition condition. 

In the Sections Q] and [3 of the paper we develop the general QS calculus of 
nondemolition input-output quantum processes in Fock space, tensored by an initial 
Hilbert space. We introduce the QS calculus of such processes using the ^-algebraic 
Minkowski metric structure of the basic quantum processes and the simple and 
convenient notation developed in [^. The Fock representation of this structure is 
closely connected with the Lindsay-Maassen kernel calculus of [S] but is given in 
terms of the matrix elements of operators for general quantum noise in Fock space 
instead of their kernels. We define the QS integrals in the framework of the new 
noncommutative stochastic analysis in the Fock scale which is described in the first 
section. 

In the Sections 01 and 0] we give complete proofs of the results, first formulated 
in [ 7 ], for the general (non Markovian) quantum filtering from the viewpoint of QS 
calculus. The advantage of the * - matrix notation enables us to prove the main 
filtering theorem for general output process as by using the indefinite metric for 
the corresponding * - algebra of generators of these nondemolition processes. 

The Markovian nonlinear filtering problem in the framework of quantum opera¬ 
tional (non-stochastic) approach was first investigated in [3 , and the possibility of 
deriving the stochastic equations of quantum filtering within this framework was 
shown in [3 . The Markovian filtering for the quantum Gaussian case and the corre¬ 
sponding quantum Kalman linear filter, first obtained for the one-dimensional case 
in 13 El, is considered using the QS calculus approach in |10j . 

The present paper is devoted essentially to the study of the nonlinear problem, 
extending the innovation martingale methods of the classical filtering theory iniini 
to the noncommutative set up of our problem. An application of the quantum 
filtering theory to the solution of the problem of the continuous observation of 
quantum spin states is given in Section El 
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gestions during the preparing of the paper. The first part of this paper was written 
in the Physics Department of the University of Milan and the second part in Gentro 
Matematico V. Volterra of the University of Roma II, for the hospitality of which 
I am very grateful. 
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1. QS CALCULUS OF INPUT BOSE PROCESSES IN FOCK SPACE 

Let US denote hy T = T{£) the state space of the one-dimensional Bose-noise, 
that is the Fock space over the Hilbert space S = L^(K+) of square-integrable 
complex functions t i—s- tf{t) on the real half-line R’*'. One should consider T as 
the Hilbert space r(£) = L^(H(]R+)) of the square-integrable functions t ^{j) 
of T = (ti,.. . ,t„) with ti G K+, ti < ■ ■ ■ < n = 0, 1 , 2, ... and scalar product 
< >= S ^ij)*x{T)dT, 


[ ‘P{T)*x{T)dT = '^ f f (p{ti,...,tn)x{tl,---,tn)dti...dtr, 

^_Q 1 ^ 


where the integral is taken over the set H(]R+) of all finite chains r on K+ with 
respect to the natural Lebesgue measure dr = dti... dt„ for every n = |t| = 
0,1,.... Following [3 we shall identify the chains r = (ti,..., with the finite 
subsets {ti,... ,tn} C M'*', so that the empty chain (n = 0) is identified with the 
empty subset r = 0 having dr = 1 and r = t (n = 1) is identified with the one- 
point subset {t} having dr = dt. We shall also denote the normalized vacuum 
function as the Kronecker (5-function: 5%{t) = 1, ii t = 0; <50(r) = 0, if r ^ 0, 
and consider the Fock spaces Tg = r(fs), = r(ff) over orthogonal subspaces 
£s = = 0 : t < s},£^ = {<^{r) = 0 : r ^ as the function Hilbert spaces 

L^(Os), LF‘{£11) on the subsets fig = {r c]s,oo[}, = {r c]t, s]} of the chains 

T > s, s > T > t correspondingly. 

Note that for any t > s a chain r G il can be represented as the triple r = 
(r*, T^,Ts) of the subchains rs = {t^ G r : ti > s}, rf = {ti G r : s > ti > t}, r* = 
{ti G r : ti < t} so that the direct product representation H x x Hg 

holds and, hence, the tensor representation T = T* ® ® Tg with iF* = I? (H‘), 

= {r G H : r < t}. 

The basic processes for QS calculus in Fock space T are the annihilation H_, 
creation and quantum number N processes, represented for all t > 0 by the 
unbounded operators 

(H_(t)v5)(r) = f y)(rUs)ds, 

Jo 

iA+{t)ip){T) = ^x\s)ipiT\s), 

sGt 


with the common dense domain = {(^j G iF : J |r‘||(i9(T)pdt < oo} , and 

(1.1) {N{t)(p){T) = |r*|<p(T), 

where x‘(s) = 1, if s < t, x*(s) = 0, if s > t, |r*| = chain 

T U s is defined almost everywhere as (ti,... ti, s, ti+i,..., t„), if ti < s < ti+i, and 
t\s = (ti,..., ti_i, ti+i,... tn), if s = t, t\s = t, if s ti for all i. Note that the 
processes H_, A'^ and N are non-commuting, but commuting with increments: 

(1.2) [H_(t), ^■'■(t')] = tAt'I, where t A t'= min(t, t'), 

[H_(t),iV(t')] = Gl_(tAt'),[tV(t),H+(t)]=H+(tAt'), 

the processes A- and are mutually adjoint: A^(t) = A-(t)* = A'^(t), and N 
is selfadjoint: N* = N. 
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Let US introduce the notations [S] 

(1.3) Al:it)=tl, A°_{t) = A_{t), At = A+{t),Al{t) = Nit), 


where / is the identity operator in J-, thus defining a 3 x 3 matrix-valued QS 
process A = (A^), indexed by /r, p G {—, o, +} with A)) = 0, if M = + or p = —. We 
shall consider the process A defined as a linear operator-valued function A(c, t) = 
tr{c A(t)} in terms of a 3 x 3 - matrix c = (c))), 

(1.4) A(c, t) = /c_|_t -|- A_(co , t) + A“''(c^, t) + N(c°, t), 

where A_(c“) = c“A_, A+(c^) = c^A+, A(c°) = c°N, writing the matrix trace 
as tr{cA} = c^A)) by the tensor notation of the sum ^ c^A)). The matrices c with 
cO = 0 for /i = -|- or p = — form a complex Lie x-algebra with respect to the matrix 
commutator and the involution 

(1.5) 

/ 0 c- c; \ / 0 c- c;* \ / 0 0 1 \ 

c = 0 c° c+ 1 -^ c* = 0 c°* Co * = gcfg, g = 0 1 0 , 

\000/ \00 0/ \100/ 

where c))* = c~^ = cf and g^ = g = g~^ is the indefinite metric matrix, defining 
a pseudo-scalar product in : 

(x|z) = x'^z~ + x°z° + x~z^ = x*z , 

X* = {x+,Xo,x-) = xfg is the row, conjugate to the column x = (x^) G C^. 

Now we can consider a multi-dimensional Bose noise, when £ is a Hilbert space 
L^(R+ ^ C™) of vector-functions (p{t) = (<P'^)(t) = (p°(t),j = 1,..., m with 


< <f\(p >= 



It is enough to regard c^ as a m-row with components Cj G C, c^j_ as a m-column 
with components gC, and c° as a m x m-matrix with elements c^ G C. The 
following theorems are valid also for the general situation T = r(L^(R+ ^ /C)), 
if the indices p,, p take values in the set —, J,-I-, where the one-point index value 
/r, p = 0 is split into to = | J| points j G J of an index set J for a basis in a Hilbert 
space K, with the infinite cardinality | J| = dim/C. 


Proposition 1. The basic QS process A{c), defined by iH.gj) . gives for eacht 

an operator representation of the complex Lie*-algehra of matrices i 1 .5\} : A(c, t)* = 
A(c*,t), 

(1.6) [A(c*, t), A(c, t')] = A([c*, c],t/\t') . 

The multiplication table ^ for Ito differentials (iA_,dA+,dA, and Idt can be writ¬ 
ten in terms o/A(c, dt) = ctdAf{t) as 

(1.7) A(c*, dt)A(c, dt') = A(c*c, dt n dt'), 

where dt P| dt' = 0 for t t', A(-, 0) = 0, and dt P| dt' = dt for t = t'. 

Proof. Taking into account, that A1 = A+ and N* = N, one obtains 

(1.8) A(c, t)* = Icf*t + A_(c^*, t) + A+(c-*, t) + A(cr, t) 

The comparing of with gives the *-property A(c)* = A(c*) of the map 
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The Lie representation property follows directly from the canonical commutation 
relations 


[A_ib-,t),A+idl,t)]=tb-d°+, [N{bl),Nidt)]=N{[bl,dl]) , 
[N{bt),A+{dl)] = A+{bld%), [A_{b:),N{dl)] = A_{b-d:) , 

which give [A(b),yl(d)] = A([b,d]), where we take into account that 

(b d)- = b-di, (b d)- = b-di (b d)^ = (b d): = bid: 
for matrices b, d of the form dni). 

Applying it to b = c*,d = c and taking into account the commutativity of 
A{c*,t) with increment A(c,t') — A{c,t), one obtains 11 .HU . In the same way one 
obtains 113 from the Hudson - Parthasarathy multiplication table 

dA_{b-)dA+{d%) = Idt{b-dl), dA_{b-) dNidl) = dA_(&-<) 

( 1 . 9 ) dN{b:)dA+{di) = dA+{b:d%), dN{b:)dN{d:) = dNibidi), 

for dA^{t) = A^{t + dt) — A^{t),h = c*,d = c. Due to complex linearity of the 
map c A(c) the formulas 11 .bll . 11.711 can be always extended to arbitrary b, d 
by polarization formula 

3 

A(bd) = ^A((b* + rd)*(b* + rd))/4i” , i = 

n—0 
3 

A(b)A(d) = ^ A(b*+ i”d)*A(b*+ i”d)/4i”. 

n—0 

Hence, 11.611 is equivalent to Ijl.SII and Ijl.TII to Ijl.bll . ■ 

Let us now define a QS integral with respect to the basic process A for a ma¬ 
trix quantum process C(t) = {Cj:){t), G {—,J,-I-} in IF. Assuming that 
the operator-valued functions t C^{t) are weakly measurable and adapted: 
C{t) = C* ®It, where C* are the operators in T* for all n G {—, J} and v G {J, 
one can define in the case of finite J the QS-integral 

f A{CAs) ■■= r^c/fdAC^ fc'idA'; 

Jo Jo Jo 


as the sum of the Lebesgue operator-valued integral f C_i_ (s)ds and the Ito integrals 
/ C-dAi, f CidA+, f C^dA^f in the sense 

In the general case £ = L^(]R+ ^ /C) we shall regard the QS-integral f Cj^dAI 
as a continuous operator 1F+ —> T- on the projective limit 1F+ = C\n>i^id) info 
T- = ClnKi^id) of Hilbert spaces G{C) C IF c C > 1 > Cj with respect to 
the scalar products 


M\v) 



((p(r)|(p(r)) + 

(F('r)iF('r))- 


dr , 


(P{t) G£(t) ,?7 > 1 
(fir) G £'{t) ,t] <1 . 


Here ((^|(/3)+(r) > ||(/3|p > {(p\(fi)-{T) are the square-norms in the Hilbert tensor 
products £{t) = ‘S>teT£{t), 'S'(r) = 'S>teT£'{t) of Hilbert spaces £{t) C 1C C 

£'{t), forming a Gelfand triple for each t G with respect to the scalar product 
IIfP = (fIf) in a Hilbert space /C (or simply £{t) = JC = £'{t), if /C = C™). 
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We shall say that a weakly measurable function 1 C(t) is locally QS-integrable 
if its components C/f, ^ G {—, o}, v G {o, +} are locally LP-integrable as operator¬ 
valued functions 


c+{t):g+ 

ciit):g+ ^g_(g,s'it), 
c-it):g+(^£{t)^g-, 
C°{t) :g+<»£it)^g_<^£'{t), 


iic+(-)ii«’:<oo (p=i) 
||C;(.)||f?<oo (p = 2) 

b = 2) 

l|C°(-)ll|r <o« iP = ^) 


Here the norms are defined for any t > 0, ^ g] 0, 1[ and a sufficiently large C > 1 by 


= f\\C^{s)\\lds, ||C:||«;r = ess.<,sup||C:(s)||«, 

||c;||« = sup{||c;(^||(0/||^||(C)}, ||C:||« = sup{||C„V||(C)/||v^°||(C)}, 

<p 

where 

e c?(C), ||<^f (C) = <<P|<^>(C), 
<P°GS(C)®f(s), ||¥^“||(C) = <<p“|<p°>(C) 


and ||C° 


I?;? = r: 


II Co- 


\f" = ^o 


1^ are the norms 


lie'll!'? = ( [\\\c{s)\\irdsy/^ 

Jo 

of the operators 

ci ^: 0(C) ^ 0(C) ® (e^V) (s) = (s)^, s < t 

c-t ■■ 0(C) ® ^ 0(C), = [ C-is)ip°is)ds 

Jo 

in the Hilbert spaces £* = © f* £(s}ds, £'* = © /p £'(s)ds. 

The following theorem shows the continuity of the QS-integral of an integrable 
C, defined on 0+ even for nonadapted e^(t) by the formula 

^(C, ds)!/?) (r) = {C^{s)^ + C~{s)^°) {T)ds 

S<t 

(1.10) + £(C';(s)^ + C“(s)^:)(r/s) , 

sGt 

where ^p^ G 0^ © £ (t) is defined almost everywhere as the tensor-function (r) = 
lp{t U t). 


Theorem 1. Suppose that C{t) is a locally QS-integrable function i.e. for any 
C < 1, t > 0 there exists C > 1, such that 


C+||!;?<oo, ||C';il!;?<oo, lie. 


«.2 


' 'Ic’j < 




|^.oo 

k,t 


< oo. 


Then the QS-integral is defined as a continuous operator ^Q(C) : 0^ 

with the estimate 

V- 


( 1 . 11 ) 


H(C(s),ds) 


<r- 


Tii!’j 


+ 


1 


c: 


Icj- 


||C“ 


i!:? 


+ -iic^ 


|?,oo 

\c,t 
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for the norms ||*o(C)|l^+ = sup^{||*^(C')</j||(7?_)/||(^||(?7+)}, where r]_ < f-e, 17+ > 
^ + c and 0 < £ < f. Moreover, the adjoint integral 

A{C,As)*ip\x >=< </5| f A(C,ds)x >, G 

Jo 



is also densely defined on Q'^ C G- as the QS-integral Jp Al(C*, ds), and the func¬ 
tion C*(t) = gC{t)*g, 


c*{t)f = c^ity, c*{t)i = Cfity, cyt)f = cuty, cyty = cm* 

is locally QS-integrable with = ||C'^(i)*||^ < 00 for almost all t. 

Proof. . In order to show the continuity of the integral n~Tnii in the projective 
topology of n^>i G{C), one should prove that 

II / ^(C(s),ds)v5||(?7-) < c||79||(?7+),||79||(r?) = (< V3|79> ( 77 ))^/^ 

Jo 

for any (f G G{ri^), 'q_ < f and a 77+ > c > 0 . Due to the definition 


A{C, ds)(p 


< 


cfdAmw 


CfdAtg^W 


+11 / c:div>ii 

^0 


C+ds+ll, 


where 


<t 

CydA°_g,= fcy{smm, { /‘c';dA+^) (r)=^(C'°(s)v7)(r\s), 

sGt 


s<.t 


and 


The first two integrals in (trrnii can be easily estimated 


c:dK^)ir) = j2ic:is)mir\s). 


as 


C,(pds 


CydAfg, 


iv-) < f\\cyisMms< f\\cy{s)\\ldsMiO 

Jo Jo 


iicTii?::ii+ii(c), 


ir,_) = ||C„>°||(0<||C-J|«||^°||(C) = 


1/2 


where we took into account that 

ll+“f(C)= J I&'MrUt)rdrdt = J |r|Cl"l-i||v7(r)fdr= 
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In order to estimate the integrals of Ci and C° let us find 


CldAt^ 


iV-) = 


s<t 


^(C'+(s)(p)(r/s) 


sGt 


? 7 _'dr = 


V-fff < [C'+(si)(/5°J('r)|[C°(s2)<p°J(r) > Tylr'drdsidsa + 
Jo Jo Jci 

rt 




by Schwarz inequality. In the same way we get 


c:dN:^ 


iv-) = 


s<.t 


'^{C°{s)ip°{s)){t/s) r^lr'dr = 


sGr 


e f f [ < [C:(siX:j(r)|[C:(s2X:j(r) >,7L"'drdsids2 + 

Jo Jo Jn 

r]_ Jj\[C°{s)ip°]{T)\\^r]'l^dTds < r]_ + \\C°(p°\\'^{r]_) 


where (p°s°^s2i'^) = U si U S2). 

Taking into account that for any £>0, ^ = ?7 + e 

(??) < ^ {M^iv + £)- llv^f (??)) < ^ll^/^f (^) , 

one can find that (l + ?7 < f I|V5||?(C). 

(i + v^y\ci^\\nv) < f (r:il!;?)'ll‘^f(C), 

(i + ^^)l|CoVII?(^) < f (llc'°llf;r)'^ll^f(C), 

if £ < ^. Hence, due to \\(p{r]^) > Hv^lKC + £) > llv^lKC) for ??'''> C + we obtain 
for ? 7 _ < ry = ,^ — £, ^ < I 


£l(C, ds)</5 




I'C:* 


^(||C-||«;? + ||c;||«;?)+i||c:||«; 


if ||‘^||(??''') <1, 0 < £ < ^, what is equivalent to (II.fill . 

Due to the duality ^(C)* = 0(C~^) of J^iC) and Q{l/C,) the QS-matrix process 
C*(t) is also locally QS-integrable, and there exists the adjoint integral f* A(C, ds)*, 
defined as in inni) by C*: 


<(p\ [ A{C,ds)x>= [ < <p|C'+(s)x + C'o (s)x°(s) > ds 
Jo Jo 

< ^:\C° 4 s)x + c:{s)x: >ds= f < C+(s)> + c;(s)>:|x > ds + 

<CT(s)> + C°(s)>:|x°>ds=< [ AiC\ds)ip\x> . 

Jo 

Obviously, II/+H(C*,ds)||J/]j+ = ||/+H(C, ds)||^; 
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Corollary 1. IfC{t) are the simple measurable adapted functions, then the defini¬ 
tion coincides with the QS-integral, given by integral Ito’s sums with respect 

to the processes m- Moreover, the QS-integral is a limit of such integral 

sums in the inductive operator topology, defined by the norms EH, if locally QS- 
integrable matrix-process C can be uniformly approximated by a sequence of simple 
operator-valued processes with respect to the defined - norms on ]0, i]. 


2. QS CALCULUS OF OUTPUT NONDEMOLITION PROCESSES 

Let US consider an initial Hilbert space = t) with identity operator 1, = 

1) 0 and denote by = t) 0 CJ* and by /* = 1 0 I* the corresponding multipliers 
of the Hilbert space H = T-fi ® Qt and identity operator I = 0 I*. Let us 

identify the basic QS process A = with the process A = 1 0 A, in Ti: 

Af{t) = tl, Afi = 1 0 Afi, A^ = 1 0 A^, Al = 1 0 Nfi A QS matrix process 

C = ^vif) acting in Ti. is called adapted, if C(t) = C* 0 It, for any t, 

where C* is a matrix of operators in Ti*. We define the QS integral of an adapted 
QS matrix process C as in n~Tnii by the sum of integrals 


(2.1) ^ A (c, ds) = ^ (c+ds + CfdAf + CldAt + C^diV^) 


which exists as an adapted process with the QS differential A(C, dt) = Clf{t)dAf{t) 
for weakly measurable, locally integrable functions t i-^ Clf{t), called below QS 
integrable processes. 

Now let us consider an adapted process X{f), defined by the QS differential 
equation 


dX{t) = (^Flfit) - X{t)5>f) dA^^{t) , A(0) =x0l, 

having the solution X{t) = X 0 I Cf;dA'f, iffC = F — A0d satisfies the 
conditions for the existence of the integral 112.1 II . where X 0 8 = (^X6^^. We 

shall define the elements Fjf of matrix-operators F(t) also for p, = — = v and for 
/r = + = u by FZ = X = and assume that Fjf = d, ii pi > v under the order 
- < o < +. 


Proposition 2. If the QS process X satisfies the QS differential equation (11), 
then the process {X*X)(t) = X(t)*X{t) satisfies the equation 

(2.2) d(A*A) = (f*F-A*A 0<5)''dA(;, {X*X){0) =x*x 0 I. 

This QS Ito formula establishes an ^-algebra isomorphism from the QS differen¬ 
tiable processes X into the algebra of matrices of operator processes Fjf defined 
above. In particular, X is formally normal (selfadjoint, unitary) iff [F,F*] = 
0, ^F* = F,F* = F“^^ with repsect to the -k-operation F*(<) = gF(t)*g, and X is 
partially isometric (isometric, orthoprojection), iffYY*¥ = F (F*F = j0d,F*F = 

f;. 
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Proof. Taking into account that 

dX = d+dt + C- dl+ + C° dN°, 
dX* = d+*dt + Cl* dA°_ + d-*dAt + d°* diV°, 

and using the QS ltd formula P, defining the product {X*X){t) = {x*x) 0 J + 
/p d(X*X) by the QS differential 

d{X*X) = dX*X + X*dX + dX*dX = 

= (C*IX + X*Cl + dAl = 

({C + X®8Y{C + X®8)-X*X®5Y dAl, 

we obtain the equation 12.211 with F = X (g) ^ + C. Due to the linearity of 12.211 
with respect to the pairs (F,X) and (F*,X*), it can be extended to 

(2.3) d{X*X') = (f*F' - X*X' ^8y dAl 
by the polarization formula 

3 

X*X' = XI (^ + ^ 

n—0 

Hence, the formula E2) is equivalent to QS Hudson - Parthasarathy Ito formula 
2] and * - property 

dx* (t) = (f* {t)-x* (t) (g) <5) dH);, X* ( 0 ) = x*® I, 

which follows from it for F' = I ^ 8, corresponding to X' = I. So the map X F 
is a homomorphism with respect to the associative operator algebra structure of X 
and F with the appropriate involutions. Furthermore it is an injection, because if 
F = 0, then X = 0, as FT = X = F+. 

Conversely, if X = 0, then F^dAl = 0 for all t, but it implies F^f = 0 due to 
the independence of stochastic integrators m 

Now let us consider an adapted selfadjoint QS process Y, satisfying a QS equation 

(2.4) dY{t) = {Z*GZ-Y^8rAt)dAl{t) , Y{0)=y^I, 

where G* = G is a * - selfadjoint matrix adapted QS process with GZ = Y = 
cl, Gl = 0 for y>v, and Z = [Z^) is a * - isometric or * - unitary matrix adapted 
process: Z*Z — I ® 8 (Z* = Z“^) such, that GZZ* = G = ZZ*G (otherwise G 
should be replaced by ZZ*GZZ*). 

We shall demand that ZZ = I = Zl, Z^ = 0, if fj, > v, and Z^ y. Y + or 
V Y ~ satisfy the conditions for the existance of QS isometric (unitary) evolution 
U{t) : defined by the QS equation (??) with X = U, Fjf = UZ^x = 1: 

(2.5) dU{t) = U{t) {Zlit) - 181) , C/(0) = J. 

Sufficient conditions for this are the conditions of local integrability of the weakly 
measurable processes Z^ in the sense of the L^’-norms [H]: 

||Z;||« < 00 , ||Z°||f) < 00 , ||Zo-||f) < (X 3 , \\Z°\\[°°'> < 00 . 
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Let US call the process Y an output process, if Y is nondemolition with respect 
to the QS process Z, generating the evolution (1^ . By this we mean the commu¬ 
tativity condition 

(2.6) [r(t),X(s)] = 0, Vt<s 

with respect to the all QS processes X{t) = /i, u G {—, J, -|-} (the conditions 

are nontrivial for fJ- + and v ^ 

Theorem 2. The process Y is defined by as an adapted selfadjoint QS process 

^ff 

(2.7) UY = YU, Y{t)=y^I+ [ D^dA-;,, 

Jo 

where D is an adapted QS integrable matrix process satisfying the conditions D^U = 
UD^, = GO — Y5^. The process Y is an output QS process, iff 

(2.8) U{s)Y{t) = Y{t)Uis) , \/t<s, 

which is equivalent to the condition \Y(t), Z!f;{s)]U{s) = 0 for s < t,Zf;U = UZf;. 
The output process Y satisfies the nondemolition condition with respect to an 
adapted QS process X, defined by 

(2.9) UX = XU, X{t) =x^I+ [ GOdl)), 

Jo 

iff [y(t),X(s)]{7(s) = 0 for all t < s. The last is equivalent to the commutativity 
conditions [Y{t), Ff{s)]U{s) = 0,Vt < s, 

(2.10) [y,^] = 0, [D,F]C/ = 0, 

where F = C -|- X 0 (the conditions are nontrivial for y + and v —). 

Proof. We obtain 12.4|l with GO = U*G^U from for U, satisfying H2.5II simply 
by applying to F = U*YU the QS Ito formula (12.411 : 

d{U*YU) = (z*(C/* (g) d)G{U (g S)Z - U*YU (g dl)). 

Conversely, we obtain 

d{UYU*) = {{U (g <5)ZZ*GZZ*(G* (g ^) - UYU* (g 5)(;, dT)), 

so the QS-process Y = UYU* obviously satisfying the condition YU = UY, is 
defined as QS integral in 112.711 with GfU* due to the assumption ZZ*GZZ* = G 
for a weakly measurable, locally L^’-integrable G. 

If the processes Y and X = Zjf satisfy the commutativity conditions ESP, then 
the isometry 

(2.11) U{t,s)=T+J^ U{t,r)[zfir)-Td(;)dA;{r) 

commutes with Y{t), as it can be easily proved by induction with respect to n = 
1,2,... for the corresponding QS Ito integral sums 

n— 1 

( 2 . 12 ) Ur.it, s)=T+Y, Urit, U) [zfiU) - Ts^) (Ic (t.+ i) - A'dtr)) , 

i=0 
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where ti = t + i{s — t)/n,Uo{t,t) = I. Hence, taking into account that U{s) = 
U{t)U{t, s), we obtain 12.411 : 

U{s)Y{t) = U{t)Y{t)U{t,s) = Y{t)U{s). 


Conversely, multiplying from the left side hand by U{t)*, we obtain the com¬ 
mutativity condition for Y{t) and U{t, s), which is equivalent H2.fill for X = Zj^ due 
to the approximation of (ITTlTl and adaptedness of Y. The condition (EH) in 

the terms of Zj^U = UZ^ can be written as: 


U(s) [Y{t),Zf^{s)] 


Y(t),Zl^(s) 


U{s) = 0, Vt < s. 


In the same way the non-demolition condition for an output process Y with respect 
to a QS process X can be written as [F(t), X(s)]17(s) = 0 in terms of XU = UX. 

Representing X in the case (EH in the form of (im as the solution of the QS - 
equation 

(2.13) dX={Z^FZ-XZ)S)'^dA''^, X{0)=xZ)I, 

and taking into account that due to (EH 

(2.14) d(yx) = (Z*GFZ-yX(g)d)OdH);, 

one can easily obtain, that [F, X] = 0, iff [y, x] = 0 and [G, F] = 0 due to ZZ*F = 
F = FZZ*. In order to satisfy the condition [F(t),X(s)] = 0 for all t < s, it should 
be completed by [Y{t), FU‘{s)] = 0 at least for y -I- or p — and alH < s due to 
the QS integral representation 

X{s) = X{t) + J^ (Z*FZ - X 0 ^)0dA;) for s>t 

and commutativity of Y(t) with Z(s) at s > t. So [D, F] = [G, F] — [F (g) d, F] = 0, 
what gives the necessary and sufficient nondemolition conditions, which can be 
written in the terms of F, D,F as (imu by multiplication on the right by the 
corresponding U. ■ 


Corollary 2. The process X is an evolute transformation X = U*{x®I)U of 
an initial operator x G B{t)) with respect to a QS unitary process U, described by 
the QS equation iff it satisfies the QS equation 12.1 A) with F = X ® S. The 
process Y is an output process with respect to the QS Markovian evolution defined 
on the von Neumann algebra A = B(i)) by the transformation Zf; = U* (zjf ® I)U 


of the initial QS generators zlf, acting in (), iff 


Y{t),zfZ)I 


0 for all t and 


pL,v. The output process Y is nondemolition with respect to X = U*{x ® I)U for 
arbitrary x G H(l)), iffY = 1® B, where B is an adapted process in Fock space Q. 


Indeed, d X{t) =x ® I is a time independent adapted process, then it satisfies 
the QS equation dX = CfdAf corresponding to = 0 = Fjf — z®16^. Hence, the 
process Z = U*ZU satisfies the equation 12.1311 with Fjf = U*{x ® I)U6^ = XSf. 
The output condition Y(t), Zf{s) U{s) = 0, for Z(s) =z®I and unitary 


U, means 




= 0 for all t, p,, p; moreover the nondemolition condition 

[F(t),X(s)]17(s) = 0 for X(s) = z® I with arbitrary x G B{i)) is possible only if 
Y =1®B. 
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Note that an output QS process Y = U*YU is defined as the sum oiy® I and 
a QS - integral 

[ A{'D,ds)=U{t)* [ A{f>,ds)U{t) 

Jo Jo 

with the QS differential A(D,dt) = (Z*DZ)(^(t)dAlJ^(t). In the case of commuting 
matrix elements D" = U*D'^U and = U*ZIU, as it happens for Zl{t) = 

= 1 (g) this integral can be defined as the QS Ito integral 

f A{B, ds) = f {D^ds + D:dA°_ + D°+dAt + D^dN:) = [ D^.dA’; 

Jo Jo Jo 

with respect to output annihilation creation A'^ and quantum number N° 
processes 

Kit) = 


A°_{t) = ^ (zidAt + Z“ds) = Atity 


zfz:)idNt + {zfziydAt + {zy:z:)kdAy + {zy^zyds 

as the unitary transformation A^ = U* A^U of the input canonical processes 

3. QS NONLINEAR NONDEMOLITION FILTERING 

Let US consider a selfadjoint family Y = {Yi) of commuting output processes 
Yi,i = I,..., n, defined by 

rt 


(3.1) 


Yyt) = y,^i+ / (z*D,z);;dA;:, 


which are nondemolition with respect a QS process X: 

(3.2) [Yyt), Z(((s)] = 0, [Yyt), X(s)] = 0 Vt < s 

As follows from llTTnil for X = yk, F = + Yfe g) the family Y satisfies the 

selfnondemolition condition 


(3.3) 


[yi(<), Yfc(s)] = 0, ys,t;i,k, 


iff 


Yi(t), D^{s)k U(s) = 0 for all t < s, and 


(3.4) 


[yz,yk] = 0 , 


Di, Ds 


U = 0, yi,k. 


Let us denote by At = {F/ : i = 1,..., n}' the reduced algebra of bounded opera¬ 
tors in H, corresponding to the measurements of the process Y* = {Y(s) : s < t} up 
to a time t, defined as the commutant of all Y* = U*Y*U, and hy O = {i/i,... ?/„}' 
the initial algebra, defining Aq = O Z> The nonincreasing family {At) is the 

family of maximal von Neumann subalgebras As C At C Aq, s > t > 0 of the ini¬ 
tial reduced algebra Ao, with respect to which F is a nondemolition commutative 
vector process in the sense of the definition Yi{t) S A) Vf (or Yi{t) is affiliated 
to A() of a nondemolition QS process given in The Abelian algebra B* = A( 
generated by F‘, with = O' ® I, generated by F° = y g) I, forms the center 
B* = At f] A't of At, hence At is a decomposable algebra, having the conditional 
expectations with respect to B* = At for any normal initial state on Ao 12 At ■ 
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As follows from the next theorem, the nondemolition principle is not only suffi¬ 
cient, but also necessary for the existence of compatible conditional expectation on 
At with respect to Bt for an arbitrary initial state vector 

We shall explicitly construct the conditional expectation not only for bounded 
X G At, but also for X affiliated to At■ An operator B is said to be defined almost 
everywhere with respect to the pair if it is densely defined in the support 

subspace JC* = P^H, where P* = inf{P = P*P G B* : P^ = C}, U\\ = 1. 

Theorem 3. Let B^ G1 be CL 'VOTL Lleu'iH/QiTL'n subclyebT'CL oti cl Lhlbect space BL- 
Then a conditional expectation e* as a positive projection onto B*, satisfying the 
compatibility condition < ^|et(X)^ >=< >, for all X G A exists on At for 

an arbitrary ^ GBL, iff B*" commutes with At : B* C A'f In this case the algebra At 
can be extended to the commutant of B*, such that for any operator X, commuting 
on the domain At^ with B* = A'^, the expectation etiX) is given on At^ by 

(3.5) et{X)Ai = AEtXi, VA G A, 

where Et G At is the orthoprojector on B*^. The formula uniquely defines 

et{X) as an operator et{X)P* affiliated to B^ on AC even for unbounded X. 

Proof. Let us suppose that [X,B] 7^ 0 for an Ai G A and B G B*, and that 
et ■ At ^ B* is defined as a positive projection, compatible with ^ G TL, for which 
< B]f >7^ 0 . Then, due to the modularity property 

et{XB) = et{X)B , et{BX) = Bet{X) , 

where X G At, B G B*", we would have < C|[et(a;),P]C >=< CI[-^A]C > 7 ^ Oj what 
would be possible only if B* would be non-Abelian. But for non-Abelian B*" the 
conditional expectation does not exist for all vectors C G TC, as can be easily shown 
for a factor S* Cl. Indeed, in this case such a vector f has to be of the form 
Co ® Cl) and e(A (g) B) =< Col^^Co > h C) B, where C, if A 7 ^ Bf 

G Hi = yB*C, corresponding to the decomposition H = Ho ® Hi. So, it is 
necessary that B* C A • 

Let us define e* for such an Abelian algebra B*' by (TTKll with A't = B* and a fixed 
f G H. The orthoprojector Et commutes with A't due to the invariance of £t = A(C 
with respect to the action of the algebra A- Hence the operator EtXEt commutes 
with A'tEp. 

EtXEtBEt = EtXBEt = EtBXEt = BEtXEt = EtBEtXEt , 

if the operator X commutes with the all B G A'^. But this means that EtXEt is 
affiliated with the reduced von Neumann algebra EtAtEt on £t, coinciding with its 
commutant on £t because the induced Abelian algebra das the cyclic vec¬ 
tor C in £t. The commutativity of EtAtEt = AfEt helps to establish the correctness 
of the definition 13.511 of the linear operator et{X) on AC AC = 0 => et{X)A^ = 0 . 
Indeed, ||et(a;)AC|| = 

WAEtXfW = \\{EtA*AEty^^EtXES\ = \\EtXEt{EtA*AEty/^f\\ 

because Etf = C and {EtA* AEtf^'^ ^ = 0, if AC = 0. 

The operator et{X) : AC ^ AEtXf having the range AtEtXf C Kf = AtC, 
commutes with arbitrary A G At due to the definition 13.511 . so P*etiX) is affili¬ 
ated to P^’A'tP*, coinciding with A(P* because P* G A!tf\At = Aj, if P* is the 
orthoprojector on Kf. 
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The map X ^ f-t{X) satisfies the unital property = AEt^ = A^ due to 

^ G St, and the modularity property 

et{XB)A^ = AEtBX^ = ABEtX^ = BAEtX^ = Bet{X)A^ 

for all A G At and B G A[, and, hence, maps the algebra At on the subalgebra 
A[ C At, represented on KA. 

Now let us prove the uniqueness of the representation Hd.5ll of conditional expec¬ 
tation €t as a map onto factor subalgebra A[/AfPf = A'^P^, where P* = I—Pi G A'l 
is the support of ^ which is the orthoprojector on At^ = /C*. Due to the commu¬ 
tativity of e{X) with At we have et{X)A^ = Aet{X)^ for A G At- So we have to 
prove, that et{X)^ = EtX^. But etlx)^ G A^^, because et{X) G A^ for X G At; 
hence we should prove, that < i3^|et(X)^ >=< B^\EtX^ > for all B G A'^, which 
is a consequence of modularity and compatibility conditions: 

< B^\etiX)C >=< ^\etiB*X)^ >=< ^\B*XC >=< B^\X^ >=< B^\EtX^ > . 


Remark. Note that one should identify the factor-algebra B*P* with the space 
L°°(V‘) of essentially bounded measurable complex functions on the probability 
space V* of all observed values u* = {u(s) ■ s <t}, v{t) = {vi){t) of the commutative 
vector process stopped at t. The probability measure ^(du*) =< > 

is induced on the Borel cr-algebra of V* by the spectral resolution = J v^I{dv*). 
If = fyt Pytfj.(dv*) is the corresponding decomposition of P* G A'^, then 

(3.6) P*et{X) = < X >„* P^.fiidA) , 

Jv* 

where < X >„t=< >, and the vectors = Pv*^/\\Pv*^\\ define the 

resolution Et = |Ci,t >< Hence one should consider et(X) for an 

X G At a.s a function et{X) : ^ CPyt giving for almost all trajectories u* G V‘, 

observed up to a time t, the posterior mean values < X >„t of a QS nondemolished 
process X{t). The initial conditional expectation cq with respect to — O' 0 I 
and ^ = tp ip is given for X = x (g) 7 as e{X) (g J by 

(3.7) pe{x) = < > pvpidv) , fi{dv) = ||T(du)V'f ■ 

7 Vo 

Here the vectors = Pvtp/\\tp^tp\\, u G Vo and the decomposition e = >< 

'ip^\p{dv) for Eq = ei^ I, {pv} define the decomposition p = f‘^pyp(dv) for = 
p®I, corresponding to the orthogonal resolution y = J ul(du) on the spectrum Vo 
of the commutative family y = {%) of the initial operators V(0) = y®I. 

4 . QS CALCULUS OF A POSTERIORI EXPECTATIONS 

Now let US suppose that ^ = ip ® Sj,, the output commuting processes are 
nondemolition with respect to 

( 4 . 1 ) X{t)=x®I+ [ {Z*PZ-X®5)‘i,dA''^ 

Jo 

with X G O, and F(t) G Tt where Tt is the ^-algebra of matrix-operators F = [Flp), 
commuting with Y{s), s <t and D(t) : Tt = {FU G At : [D^, F] = 0, i = 1,..., n}. 
In the following we shall also demand that the process Y = (Y) is continuous from 
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the right in the sense A'f = ns>t what is equivalent to G A[ for all i and 

t. 

Let us denote by C the linear span of the initial operators {j/i} with the operators 
yo G Co from the commutative ideal Cq = {b G O' :< biplbip >= 0}. We also denote 
by P* the ^J-span of the operator-matrices {Di}(t) with the ideal 

Vl = {T> G At : DZ = 0 = 0+, (DZ+|DZ+) = 0} 

of the commutative ^-algebra 

Ai = {D^GA't\[B,F]=0,FGA}, 

corresponding to the kernel of the pseudoscalar product 

(Z+IZ+) =< ^|(Z*Z);C >=< eiz;gz+e > , 

where Z^ = (^Z^*, is tbe conjugate row to the column-operator Z+. Now 

we can formulate the main theorem. 

Theorem 4. Suppose that the output observed proeess 1221) is nondemolition with 
respect to a QS proeess i4A]) , and the spans C C O' and D* C are * - and * 
- algebras correspondingly. Then the posterior mean value et{X(t)) for an initial 
state vector = i/’0 Jg, !/> G t), is defined by an adapted commutative vector-process 
Kt = {nf) , k\ G A'tji = 1,... ,n, almost everywhere as an A't linear nonaticipating 
transformation of the output process Y by the stochastic ltd equation 

(4.2) detiX{t)) = etiZ*FZ)f{t)dt + KliX{t))dYfit) 

Here ei(F):,: = et{F^),eo{x 0 /) = e{x) (g) I, kAYi = ^AYi, and 

(4.3) dYfit) = dYfit) - et(Z*D,Z);(<)d<, F,(0) = y. 0 J 

are the observed martingales with respect to the filtration (ct), cind state vector 
called the innovating process for (A'l). The process Kt is defined uniquely up to the 
kernel of the correlation matrix-process 

(4.4) ft,(t) = e* (Z*D*DfeZ); (t) = CtHD^Z: + DD^iD^Zf + Dl),]it) 
by the linear algebraic equation 

(4.5) g,,it)K^ = et (Z*D*FZ); (t) - et{X{t))et (Z*D*Z); (t), 
having in the case F = X 0 corresponding to XIf) =x® I, the form 

(4.6) ft,(t)4 = et [zID:*XZI) it) + et [xDl^Zf + Zl*D-*x) (t), 

where Xf) = Xf) — efiXit)). The initial a posteriori mean value e(x) is the linear 
combination e(x) =< iplxip > +K'{x)yi, ofyi = %— < iflyiip > 1, where k = (k*) 
is defined by the equation 

(4.7) =< fi\y*xtfj >, =< fi\y*ykip >, 

with X = X— < 'iflx'ijj > 1, uniquely up to the kernel of the initial correlation matrix 

Q = (ftfc)- 

In order to prove this fundamental filtering theorem we need the following lem¬ 


mas. 
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Lemma 1. If the process X satisfies the equation then there exists such a 

martingale Mt with respect to affiliated with A'^ on Atf,, such that almost 

everywhere 

(4.8) et{X(t))=e{x)®I+ [ e«(Z*FZ);(s)ds + M*. 

Jo 


Proof. Let us define Mt on f by 

Mtf = {Et - Eo) Xm + f {Et - Es) {Z*FZ)>fdA’;^ . 

Jo 

Obviously, that Mtf satisfies the (e*,^) - martingale condition EgMtf = for 
all s < t, and 

EtX{t)f = {ex^I)^+ [ EfiZ^FZ)f{s)fds + Mt^ 

Jo 

due to Es{Z*FZ)^dAf^f = Es{Z*FZ)ffdt for ^ = -0 0 5^^. The operator Mt, 
affiliated with A't can be correctly defined almost everywhere by 

MtAf = AEtMtf = AMti, VA G At, 

as in the case of (123 for X = Mt, et{Mt) = Mf 
So, for any ^ G »4t we have 

etiX{t))Af = A ® '^0 + / Es(Z*FZ)f(s)^ds + Mtf 

= {e{x)(g)I)Af+ f es{Z*FZ)f{s)Afds + MtA^, 

Jo 

and, hence, (Oil holds on the dense linear manifold At£, of the support /C* of the 
state f on 

Lemma 2. A process Mt = (Z*DZ)0dA)) with D(t) is a martingale with 

respect to (et,^), iff et{Z*F)Z)f{t) = 0 for all t, that is almost everywhere 

(4.9) Df{t) + Df{t)et{Zl) + et{ZlYD\{t) + efiZ^ DlZlfit) = 0 

and is the zero martingale (almost everywhere), iffF>{t) G T>q, which is equivalent 
to et{Z*F)*F>Z)f{t) = 0 for all t almost everywhere, that is 

(4.10) et[{Dl + D^ZlYiDf + D^Zfm = 0 
and, hence, Df_{t) = et(,Z'YD°Zf){t). 

Proof. Due to commutativity of M(t) with At, we have to prove only that 
EtMrf = Mtf, for all r > t, iff i?t(Z*DZ)^(t) = 0 for all t. Indeed, EfiM^ — Mt)f, = 
i?((Z*DZ))i(s)^ds = 0 for all r > t iff i?((Z*DZ)^(s)^ = 0 for all t < s, which 
is equivalent to Et{Z*F>Z)f{t)f = 0 for all t due to EtEg = Et for t < s, written 
in the form (Onil for 

{Z*r>Z)f =Df+ DfZf + ZffDf + ZffDlZl . 

If Mt is a martingale, then 

et[{Mr - Mt)*{Mr - Mt)] = efiMfMr) - MfMt = 
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et(Z*D*DZ);(s)ds > 0 . 

Hence, if Mt is a zero martingale, |Mr — = 0, and et(Z*D*DZ);j](s) = 0 for 

all t > s, what is equivalent to et(Z*D*DZ);|^(f) = 0 for all t, or to 14.1 OH in view 
of 

(Z*D*DZ); = (£)+ + D°Ziy{Dl + D°Zl) . 

But this means, that (D^j. + D°Z'^)^ = 0, i.e. D(t) G Pg. Conversely if D(t) G Vq, 
i.e. if (^|(Z*D*DZ):fO = 0, then £;t(Z*D*DZ);(t)^ = 0 because 

(B^I(Z*D*DZ)+0 = ((l(Z*D*(B 0 S)DZ)+0 = 0 

for any B G Af. 

Lemma 3. Let the linear complex span of {y^} and the span of {T>i}{t) with the 
coefficients in A!^, be commutative * - and -k - algebras C and T>* up to the ideals 
Co CC and Vq C V* correspondingly. Then the locally bounded process 

(4.11) B{t) = {yo + Xlyy 0l+ f (Z*(Do + A:D,)Z)0(s)dH;(s) , 

Jo 

where^ yo G C, fp = %- < y\yiy > 1 ,DqG Pg, Df; = Df;, if {p,v) ^ (-,+), 
and Dfft) = Dfft) — et(Z*r>Z)f(t}, defined by weakly measurable locally bounded 
functions 1 1 —> A) G Ag G C, compose a weakly dense *-algebra C* in Af 

Proof. Using the QS Ito formula 12.211 . one obtains for dB = {Z*'DZ)ifdA’f with 
D(t) = Do(t) + AjD,(t)Gl?* 

d{B*B) = (Z*(D*(H 0d) + (B0 ^)*D + T>*'D)Z)^dAf 

due to the commutativity of B{t) with Zjf{t). But D*(t)D(t) G and, hence 
(T>*{B 0 S) {B 0 ^)*D + D*D)(t) G D* is an A[ linear combination of {Di(t)} 
and a G G Ug, as well as b*b G C for 6 = j/g + j/^A* G C is a linear combination of 
y^ and a b G Cq. Hence, B*B is a process of the same form as B{t), what means 
that the operators B{f) compose a *-subalgebra C* of Af The algebra 6 is a weakly 
dense in A^ because it has the same commutant Alt, as the family {U(s) : s >t}, 
and hence generates the same von Neumann algebra . 

Proof of the Theorem ^ We shall look for the martingale M, defining the 
decomposition (lOll in the Lemma^ Let us suppose, that it is a stochastic integral 
nonanticipating span 

Mt = Az*D,Z)0(s)«:dH;:(s) , tcj G Al 
Jo 

of the observable martingales 

Yyt)=Yft)-J^ e«(Z*D,,Z);(s)ds = J Z*i),ZdA , 

where Yi{t) should not be taken into account, if Dft) G I?g, as it is a zero almost 
everywhere martingale according to the Lemma|21 Due to the weak density of C* in 
Al(, proved in Lemma |21 it is sufficient to find the coefficient At) from the condition 

< f\Bitrx{t)f >=< rntyetixm > 
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for all B(t) in the form 114. 11 II . Using the QS Ito formula llOl for dB = (Z*DZ)(^dA 
where D == Dq + DjA*, one can obtain 


d < c|B*xe>=<ei(z*(S05 + D)*FZ);^>df = 

= < C|5*et(Z*FZ); + (Z*D*FZ);|$ > dt . 

On the other hand, taking into account that det(Af(t)) = €t{'Z*FZ)'^{t)dt + dMt, 
dM = (Z*DiZ)f^KMAJ;^, one can obtain 

d < ^\B*etiX)^>=<^\B*etiZ^FZ)^^>dt + 

+ < el(Z*D*Z);et(X) + (Z*D*D,Z);K*|C>dt . 

Hence, < C|Z*D*D,Z|KjC >:^=< {(Z*D*FZ); - (Z*D*Z); et(X)} ^ >, what 

is equivalent for D = Dq + A* to (14.511 and 

(4.12) < C|(Z*D*D,Z);i^jC >=< C|{(Z*D*FZ); - (Z*D*Z);et(X)K > 

due to D*Di = D*Di and arbitrariness of X] € A'^. But the left hand side of the 
last equation (mUi due to Schwarz inequality is zero: 

< ci(z*D5DoZ);^ >=< {D^zi + diu\{d:z°+ + DIU >= 0 

as < <e|(Z*D5DoZ);^ >= \\{D°Zl + = 0 for Do S On the other 

hand, taking into account, that 

< e|(Z*D5FZ);^ >=< e|(Z*Do(X O S)Z)-i >, 

as < C|(Z*D*CZ);^ >=< {D°Zl + Dl)o^\{C°Zl + Cl)C >= 0 for C = F-X^S, 
and due to Tk- - normality 


D°D°* = D°*D° D°D~* = D° 


D-Dl 


Dl*D 


O 

+ 


of Do G Dq as for a matrix-operator of the commutative matrix ^-algebra Jy, one 
can obtain 


< C|(^*D5FZ);C >=< + D-zi)*x^ > + 

+ < [Dizi + Diu\xzi^ >=< + ZI*XDTZIU >= 

= < mV + Dl*D°+Dl)oXi >=< ^|(D; + Dl*D°+Dl)let{X)^ > . 

Here is quasi-inverse conjugate matrix-operator for normal = D°* D°^ D°, 

and we used 

{DlZl + DIU = 0, {DTZ°+ + = 0 

as for -k - normal Dq G Dq. Hence, the right side of equation 114.1211 is also zero: 

< ^|(z*D5Z);e >=< ?|(z*D5Z);e*(X)c >, 

because in the same way one can obtain 

<^|(Z*DSZ);et(X)e> = <CI(D;+D-Z“)Set(X)C> + 

+ < {D°Zl + D°^U\Zlet{X)^ > = ^\{D^ + 0^0°+DD^etiX)^ > . 

This proves also the uniqueness of the solution of the equation (lOll up to the kernel 
of the correlation matrix g{t) = (t) because if ft^Ao = 0 for an A'^. - adapted 


■p: 
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vector process Ao = (Aq), then Dq = AgDi G Vq, and, hence, 

< A^^l {e(Z*D*FZ); - e{Z^T>*Z)-e{X)}^ > = 

+ <^|{(Z*D5FZ);(Z*DSZ);e(X)}^> = 0. 

In the case Fj^ = XS^ taking into account that 

e(Z*D(X 0 <5)Z); = D-e{X) + D-e^XZ^) + e[Z°^X)D\ + ) 

and e(Z*DZ);e(Z) = 

Dle{X) + DZe{X)e{Zl) + e(Z;)*e(X)7A; + )e(X), 

one can easily obtain the equation EH) from EH- 

One should look for the initial condition eg (a; O /) = e{x) O I for the equation 
E3 in the linear form e(a;) =< il)\x’4> > +yiK*, where k* should be found from 
< b'il)\zip >=< hip\e{z)ip > for all 6 = yo + Sy^A*, where yo G C and A* G C. This 
gives the initial equation EH- 

Corollary 3. If {y^} are commuting orthogonal projectors in f), and also D* = 
Di = D| are commuting -k-matrix projectors, then the conditions of the Theorem^ 
are fullfilled, and they are fullfilled also in the case iAT = 0. In particular, for the 
case D° = I Z) S°, Df = 0 = D^, Df_ = 0, corresponding to the counting output 
process Y = N, the eguation EH gives 

m = et{Zl*XZl){t)/et{Zl*Zl){f) - et{X){t) , 

if €t{Z^ Zf){t) 7 ^ 0. In the other case Df = 1 = iA^, D° = 0, Df = 0, corre¬ 
sponding to the output coordinate observation Y = Q, one obtains 

(4.13) K* = et{XZl + Zl*X) - et{X)et{Zl + Z^*) . 


Indeed, the linear, span of commuting orthoprojectors {yt}, and also ^(-span 
of ^-projectors {Di} is a *-and *-algebra C and TA* correspondingly. In the case 
ZA-° = 0 the product D* • is in lAg as matrix with (D*Dfe)° = 0, (D*Dfc)5j_ = 0, 
(D*Dfc)- = 0. Such commutative matrices also form a commutative *-algebra 
T>* up to the ideal Vq, because G*G = 0, and, hence, (GZ+IGZ+j = 0 for a 
matrix-operator G with Gf = Q for (y, p) ^ 

5. An application of the QS filtering 

The applications of the filtering equation EH to the derivation of a posteriori 
Schrodinger equation for the coordinate observation are given in HTnisi, and for 
the counting observation are given in dEOl. 

In contrast to the usual Schrodinger equation, describing a closed quantum sys¬ 
tem without observation, these new stochastic wave equations give the dynamics of 
an open quantum system undergoing the nondemolition measurements which are 
continuous in time. Thus the continual wave packed reduction problem is solved 
by the quantum filtering method for the typical QS models of observation such as 
a quantum particle is a bubble chamber |1X| (diffusive observation) and an atom 
radiating the photons m (counting observation). Here we consider another ex¬ 
ample of the quantum nonlinear filtering - the QS spin localization, describing the 
continuous collapse of the vector polarization p = {pi,P 2 ,P 3 ) for the spin ^ of an 
electron under a continuous nondemolition measurement in a magnetic field. The 
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polarization p(t) at the time instant t > 0 is given by the conditional expectations 

(E31) 

(5.1) pjit) = et{X,{t)) , X,{t) = U*{t){x, 0 I)U{t) , 

where Xj = aj are the Pauli matrices 



and [/(t) is a QS unitary evolution in the Hilbert space 0 T. Here is 


the Hilbert space 


of the spinors = 


i/i- 

V'+ 


(r), where ' 0 ^(r), 


r G are the wave functions of the nonrelativistic electron with the definite 
^-projections of its spin s = ^(di, 0 - 2 , (T 3 ) and probabilistic normalization 
IIV'IP /= 1 , where = |_p -|- |'0+P) and X = r(£) is the Fock 
space over the Hilbert space E = L^(]R+) 0 C". 

The initial polarization p{0) = = POj 



i/;(f)'l'o-jT/i(f=Odf , j = 1, 2, 3 , 


has the values in the unite ball K = {p G : |p1 < 1}, where |p|^ = {p,p) = p^, 
i.e. can be mixed 1 in the pure (vector) state tp G IIV’IP = 

J '(/i(r)’l’'(/i(r)dr = 1. 

Let us suppose that the evolution U {t) defines the system of Langevin equations 
(im|l of the form 

- n n 

( 5 . 2 ) dX + ii[X,H] + -Y,[[X,L,],L,])dt = iY.^X,U]dV, . 

i=i j=i 

Here X{t) = {Xi, X 2 , X 3 ){t), H{t) = 5 spin-Hamiltonian, 

corresponding to the magnetic tense u{t) = (m^, M^)(t) G R^, 

1 " 1 

Lzit) = - '^rl{t)Xj{t) = -Ri{t) 

i=i 

are spin-operators, defined by the real vectors fi{t) = (r^,rf G R^, i = 

1,... ,n, and Vj = 1 0 i = 1,... ,n are the independent standard Wiener 

processes, represented by the input operators \{A'l{t) — A^_{t)), i = \^-A in the 
Fock space J- with respect to the initial vacuum state Jg G T. The stochastic system 
of the operator equations 15.211 corresponds to the unitary Markovian evolution (1^ 
in () 0 IF, f) = with the generators 


ZP;{t) = U{tr{ZP;{t)®l)U{t) 


p = -,%,...,n 
u = 1,..., n, -I- 


defined by the spin-operators 


1 


Zfc = (5fcl , z\ = -fi = z^ , =- 


1 (I 


2 V4 


m 


, i = , 
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where fi{t) = H*) 

evolution realizes the output coordinate processes 


3 

u^{t)aj. Such the 
i=i 


Y,{t)=U{t)*W,{t)U{t)=Q,{t), i = 
satisfying the QS equations ifOll in the form 

(5.3) dFi = R,dt + dWi , dWi = i® TRAf 


of the indirect non-demolition observation of the noncommuting spin-operators 


Ri{t) = U{t)*{fi{t) 0 I)U{t) , i = l,...,n . 


The standard Wiener processes Wi, i = 1,... ,n, represented by the commuting 
operators Af (t) + Al_ (t) in iF, describe the independent errors Vi — Ri as the white 
noises Wi. They do not commute with the white noises Vi of the perturbations in 
the quantum system (15.211 : 

(5.4) ['l/,(s) , Wfc(t)] = 2W(s - t)d,J , 

due to 

[^^(s) , Wk{t)] = 2i[y4!_(s), A+{t)] = ■ 


Proposition 3. Under the given assumptions the a posteriori spin polarizations 
EIli satisfy the following system of nonlinear stochastic equations 

^ n n 

(5.5) dp+{pAu+- '^{rlp- {p, ri)f,))dt = - {p, ri)p)dV , 

i=l i=l 

where dl)(t) = dYi{t) — {p{t), ri{t))dt. 


Proof. Let us consider the nonlinear filtering equation (lOl for the spin-operators 
d — 2, 3, which are equivalent to the Pauli matrices dj, j = 1, 2, 3. We can 

use (lOl for the evaluation of the expectations because the conditions of the 
Theorem 0] are fulfilled (see the Corollary of sec. Ej). The innovating martingales 
dli = dYi — €t{ri)dt in this case are given by the differences dl) = dF^ — {p, ri)dt 
because 

3 3 

it{Rt{t)) = '^rl{t)et{Xj{t)) = '^rl{t)pj(t) . 

i=i 

Due to Pii.{t)dt = et{dYi{t)dYk{t)) = Sikdt, the coefficients Kl{Xj{t)) are given by 

= ^f^t{Xit)Ri{t) + Ri{t)X{t)) - etiX{t))et{R^{t)) = 

= n{t) - {p{t), n{t))p{t) , 

3 . 

because djfj + fjdj = 2r^l for fi=Y ''~i^j^ 

t=i 

XjRi + RiXj = U{t)*{djfi + ridj)U{t) = 2r{l . 

The vector-product p{t) A u{t) in (15.511 represents the expectations 


iet{[X{t) , i?(t)]) =iet i^[X(t) , X,{t)]u,{t) 


i=i 
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because [ct, u] = = jS A u, and 

[X{t) , H{t)] = \u{tr{[d,u{t)]®I)U{t)- 

In the same way one can obtain 

rUt)p{t) - {p{t), n{t))n{t) = {p{t) A n{t)) a n{t) 

as for the vector representation of the double commutator , Li{t)] , Li{t)], 

defining together with i[X{t) , H{t)] the products {Z*{t)X{t)Z{t))^ in (14.211 . ■ 

Now we can prove, that the continuous indirect nondemolition measurement (Oil 
of the quantum spin reduces any initial state of the electron at the limit t ^ oo, to 
the completely polarized one. This gives a kind of the stochastic ergodicity property 
of the nonlinear system of quantum filtering equation (15.511 . 

Theorem 5. Let p(0) = po G B he an arbitrary initial polarization for the non¬ 
linear quantum filtering equation \5.5il . Then this equation has a unique stochastic 
solution pit) G B, and p^lf) = {p{t) , pit)) ^ \ at t ^ oo almost surely, if 

Ht) = fo E ki(s)Pds ^ oo. 

Proof. The vector stochastic equation (Oil up to a renormalization f(t) = 
p{t)p{t) is equivalent to the linear stochastic equation 

^ u n 

(5.6) df+{fAu+-'^{r^f-{f,ri)ri))dt = p'^ndY, . 

Indeed, let p{t) be the stochastic Ito’s integral 

(5.7) p{t) = l+ / ^(/(s), ri(s))dy*(s) 

•^0 i=l 

defined by the unique solution f{t) of this ordinary linear stochastic differential 

equation with the initial nonstochastic vector /(O) = po. Then dp = ^(/,r))dli, 

2=1 

and by Ito’s formula 

d{pp) = dpp + dpdp + pdp 

we obtain the equation for / = pp iff p{t) satisfies the equation (15.511 : 

d/+ if A u + I ^ (/, ri)n))dt = dpp + dpdp+ 

^ 2=1 

u u n 

+ - (P) ri)f))dYi = ^(/, ri)pdYi + ^(r) - (p, ri)p){f, f,)dt 

2=1 2=1 2=1 

n n 

+ X! ^ X! ■ 

2=1 2=1 

This the unique solution of the nonlinear filtering equation (15.511 with p(0) = po can 
be written almost surely (p(t) 0) as plf) = f{t)/p{t), where f{t) is the solution 

of the linear equation (15.611 with /(O) = pb, and p{t) is the integral (15.711 . 

In order to prove that almost surely |p(t)| < 1, if |po| < 1, it is sufficient to show, 
that 

fit) = (fit), ft)) < p(t)^ if /(O) = Po . 
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Using the Ito’s formula we obtain 

U 

df = 2(/, d/) + (d/, d/) = 2p^(/, rl)dr, - 

U U 

- - (/’ = V + (p" - f) Y. ’ 

2^1 2^1 

where dp^ = 2pdp + (dp)^ =2pJ2if , n){dY^ + (/,f*)). Hence 

2=1 

d(/^ - P^) = A(p2 - /^)dt, 

U 

where A = ^ rf > 0, and 

2=1 

P^W - (/W. /W) = e“^d)(i _ vt . 

Thus P[t) < p‘^{t), if |p°| < 1, and f^{t) p^{t) exponentially at t ^ oo, if 

A(t) —> oo {P(t) = p^(t) , Vt, if IpoI = !)• This proves that p{t) = f{t)/p(t) —> 1 
almost surely {p{t) ^ 0) due to the positivity of p{t). ■ 

Remark. The model (15.211 of continual nondemolition measurements of noncom¬ 
muting spin-operators Ri{t), i = 1,... ,n in the quantum stochastic system (15.211 is 
unique in the Fock space !F = r{£) over the minimal Hilbert space S = L^(R+)(8)C”. 
It can not be realized in the framework of classical probability theory due to the 
noncommutativity (lOll of the quantum stochastic processes Vi{t) and Wi{t) though 
each of them can be described as the classical one separately due to the selfnonde¬ 
molition (commutativity) property \Vi{t) , 14(s)] = 0 = [Wi{t) , Wk{s)]. 

The result obtained here in a rigorous mathematical way corresponds to a rather 
intuitive physical picture of the continual spontaneous collapse of the quantum spin 
under the non-demolition observation. This proves the appropriateness of the given 
quantum stochastic setup for the theory of continuous measurements and quantum 
filtering. 
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